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I explore the form of the effective interaction in harmonic-oscillator-based effective theory (HO- 
BET) in leading-order (LO) through next-to-next-to-next-to-leading order (N^LO). As the included 
space in a ffOBET (as in the shell model) is defined by the oscillator energy, both long-distance 
(low-momentum) and short-distance (high-momentum) degrees of freedom reside in the high-energy 
excluded space. A ffOBET effective interaction is developed in which a short-range contact-gradient 
expansion, free of operator mixing and corresponding to a systematic expansion in nodal quantum 
numbers, is combined with an exact summation of the relative kinetic energy. By this means the 
very strong coupling of the included (P) and excluded (Q) spaces by the kinetic energy is removed. 
One finds a simple and rather surprising result, that the interplay of QT and QV is governed by a 
single parameter k, the ratio of an observable, the binding energy to a parameter in the effective 
theory, the oscillator energy Uuj. Once the functional dependence on k is identified, the remaining 
order-by-order subtraction of the short-range physics residing in Q becomes systematic and rapidly 
converging. Numerical calculations are used to demonstrate how well the resulting expansion re- 
produces the running of H"^^ from high scales to a typical shell-model scale of Sfiw. At N"^LO 
various global properties of H'^^^ are reproduced to a typical accuracy of 0.01%, or about 1 keV, at 
9>hhj. Channel-by-channel variations in convergence rates are similar to those found in effective field 
theory approaches. 

The state-dependence of the effective interaction has been a troubling problem in nuclear physics, 
and is embodied in the energy dependence of W^^^ {\E\) in the Bloch- Horowitz formalism. It is 
shown that almost all of this state dependence is also extracted in the procedures followed here, 
isolated in the analytic dependence of W^^^ on Thus there exists a simple, Hermitian H^^^ that 
can be use in spectral calculations. 

The existence of a systematic operator expansion for H'^^^ , depending on a series of short-range 
constants augmented by k, will be important to future efforts to determine the HOBET interaction 
directly from experiment, rather than from an underlying NN potential. 

PACS numbers: 21.30.Fe,21.45.Bc 
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I. INTRODUCTION 

In nuclear physics one often faces the problem of determining long-wavelength properties of nuclei, such as binding 
energies, radii, or responses to low-momentum probes. One approach would be to evaluate the relevant operators 
between exact nuclear wave functions obtained from solutions of the many-body Schroedinger equation. Because 
the NN potential is strong, characterized by anomalously large NN scattering lengths, and highly repulsive at very 
short distances, this task becomes exponentially more difficult as the nucleon number increases. Among available 
quasi-exact methods, the variational and Green's function Monte Carlo work of the Argonne group has perhaps set 
the standard |T], yielding accurate results throughout most of the Ip shell. 

Effective theory (ET) potentially offers an alternative, a method that limits the numerical difficulty of a calculation 
by restricting it to a finite Hilbert space (the P- or "included" -space), while correcting the bare Hamiltonian H (and 
other operators) for the effects of the Q- or "excluded"-space. Calculations using the effective Hamiltonian H'^-^f 
within P reproduce the results using H within P+Q, over the domain of overlap. That is, the effects of Q on P-space 
calculations are absorbed into P{H'^f^ — H)P. 

One interesting challenge for ET is the case of a P-space basis of harmonic oscillator (HO) Slater determinants. 
This is a special basis for nuclear physics because of center-of-mass separability: if all Slater determinants containing 
up to N oscillator quanta are retained, H'^ff will be translationally invariant (assuming H is). Such bases are also 
important because of powerful shell-model (SM) techniques that have been developed for iterative diagonalization 
and for evaluating inclusive responses. The larger P can be made, the smaller the effects of H"^^^ — H. If one could 
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fully develop harmonic-oscillator based effective theory (HOBET) , it would provide a prescription for eliminating the 
SM's many uncontrolled approximations, while retaining the model's formidable numerical apparatus. 

The long-term goal is a HOBET resembling standard effective field theories (EFTs) [HIS]- That is, for a given choice 
of P, the effective interaction would be a sum of a long-distance "bare" interaction whose form would be determined 
by chiral symmetry, augmented by some general effective interaction that accounts for the excluded Q space. That 
effective interaction would be expanded systematically and in some natural way, with the parameters governing the 
strength of successive terms determined by directly fitting to experiment. There would be no need to introduce or 
integrate out any high-momentum NN potential, an unnecessary intermediate effective theory between QCD and the 
SM scale. 

One prerequisite for such an approach is the demonstration that a systematic expansion for the HOBET effective 
interaction exists. This paper explores this issue, making use of numerically generated effective interaction matrix 
elements for the deuteron, obtained by solving the Bloch-Horowitz (BH) equation for the Argonne vis potential, an 
example of a potential with a relatively hard core (< 2 GeV) i4j. The BH H'^^^ is a Hermitian but energy-dependent 
Hamiltonian satisfying 



Here H is the bare Hamiltonian and E and \E' are the exact eigenvalue and wave function (that is, the solution of 
the Schroedinger equation in the full P + Q space). E is negative for a bound state. Because H'^-f^ depends on the 
unknown exact eigenvalue E, Eqs. ([T|) must be solved self-consistently, state by state, a task that in practice proves 
to be relatively straightforward. If this is done, the P-space eigenvalue will be the exact energy E and the P-space 
wave function ^'p will be the restriction of the exact wave function ^! to P. This implies a nontrivial normalization 
and nonorthogonality of the restricted (P-space) wave functions. If P is enlarged, new components are added to the 
existing ones, and for a sufficiently large P space, the norm approaches ones. This convergence is slow for potentials 
like Wis I with many shells being required before norms near one are achieved [51 [5] . Observables calculated with the 
restricted wave functions and the appropriate effective operators are independent of the choice of P, of course. All of 
these properties follow from physics encoded in H'^^^. 

In HOBET P and thus H'^^^ are functions of the oscillator parameter 6 and the number of included HO quanta 
hp. In this paper I study the behavior of matrix elements {a\H'^^^l3) generated for the Argonne Wig potential, as 
both b and Ap are varied. In particular, Ap is allowed to run from very high values to the "shell-model" scale of 8 
fiLo, in order to test whether the physics above a specified scale can be efficiently absorbed into the coefficients of some 
systematic expansion, e.g., one analogous to the contact-gradient expansions employed in EFTs (which are generally 
formulated in plane wave bases). There are reasons the HOBET effective interaction could prove more complicated: 

• An effective theory defined by a subset of HO Slater determinants is effectively an expansion around a typical 
momentum scale q ^ 1/6. That is, the P-space omits both long- wavelength and short- wavelength degrees of 
freedom. The former are connected with the overbinding of the HO, while the latter are due to absence in P of the 
strong, short-range NN interaction. As any systematic expansion of the effective interaction must simultaneously 
address both problems, the form of the effective interaction cannot be as simple as a contact-gradient expansion 
(which would be appropriate if the missing physics were only short-ranged) . 

• The relative importance of the missing long-wavelength and short-wavelength excitations is governed by the 
binding energy, \E\, with the former increasing as \E\ — s- 0. These long-range interactions allow nuclear states 
to de-localize, minimizing the kinetic energy. But nuclei are weakly bound - binding energies are very small 
compared to the natural scales set by the scalar and vector potentials in nuclei. One concludes that the effective 
interaction must depend delicately on \E\. 

• An effective theory is generally considered successful if it can reproduce the lowest energy excitations in P. But 
one asks for much more when one seeks to accurately represent the effective interaction, which governs all of 
the spectral properties within P. The HO appears to be an especially difficult case in which to attempt such a 
representation. The kinetic energy operator in the HO has strong off-diagonal components which raise or lower 
the nodal quantum number, and thus connect Slater determinants containing Ap quanta with those containing 
Ap ± 2. This means that P and Q are strongly coupled through low-energy excitations, a situation that is 
usually problematic for an effective theory. 

All of these problems involve the interplay, governed by |P|, of QT (delocalization) and QV (corrections for short- 
range repulsion) . The explicit energy dependence of the BH equation proves to be a great advantage in resolving the 
problems induced by this interplay, leading to a natural factorization of the long- and short-range contributions to the 



H'ff ^ H + H 




\^p) = {i-Qm. 



(1) 



3 



effective interaction, and thereby to a successful systematic representation of the effective interaction. (Conversely, 
techniques such as Lee-Suzuki ^ will intermingle these effects in a complex way and obscure the underlying simplicity 
of the effective interaction.) The result is an energy-dependent contact-gradient expansion at N'^LO that reproduces 
the entire effective interaction to an accuracy of about a few keV. The contact-gradient expansion is defined in a way 
that is appropriate to the HO, eliminating operator mixing and producing a simple dependence on nodal quantum 
numbers. The coefficients in the expansion play the role of generalized Talmi integrals. 

The long-range physics residing in Q can be isolated analytically and expressed in terms of a single parameter, 
K — ^j2\E\/Tiuo, remarkably the ratio of an observable {\E\) to a parameter one chooses in defining the ET. The 
dependence of H'^^^ on n is determined by summing QT to all orders. The resulting H'^^f is defined by k and by the 
coefficients of the short-ranged expansion. 

This same parameter governs almost all of the state dependence that enters when one seeks to describe multiple 
states. Thus it appears that there is a systematic, rapidly converging representation for H'^^^ in HOBET that 
could be used to describe a set of nuclear states. The short-range parameters in that representation are effectively 
state-independent, as the state-dependence usually attacked with techniques like Lee-Suzuki is isolated in k. 



II. LONG- AND SHORT- WAVELENGTH SEPARATIONS IN H^^^ 

In Refs. OIS] a study was done of the evolution of matrix elements {a\H'^f , for the deuteron and for ■^He/'^H, 
from the Ap oo limit, where H'^^f H, down to Ap characteristic of the shell model (SM), e.g., small P spaces 
with 4, 6, or 8 huj excitations, relative to the naive Is-shell ground state. As noted above, this definition of P in 
terms of the total quanta in HO Slater determinants maintains center-of-mass separability and thus leads to an H'^-^^ 
that is translationally invariant, just like H. Indeed, the HO basis is the only set of compact wave functions with this 
attractive property. 

But this choice leads to a more complicated ET, as P excludes both short-distance and long-distance components 
of wave functions. This problem was first explored in connection with the nonperturbative behavior of H'^f^: the 
need to simultaneously correct for the missing long- and short-distance behavior oi^ p is the reason one cannot tune 
P to make H'^^^ converge rapidly. For example, while it is possible to "pull" more of the missing short-range physics 
into P by choosing a small 6, this adjustment produces a more compact state with very large Q-space corrections to 
the kinetic energy. Conversely, one can tune b to large values to improve the description of the nuclear tail, but at the 
cost of missing even more of the short-range physics. At no value of h are both problems handled well: Fig. [T] shows 
that a poor minimum is reached at some intermediate &, with a 10?ia; "bare" calculation failing to bind the deuteron. 

The solution found to this hard-core correlation/extended state quandary is an a priori treatment of the over binding 
of the harmonic oscillator. The BH equation is rewritten in a form that allows the relative kinetic energy operator to 
be summed to all orders. (This form was introduced in the first of Refs. [B]; a detailed derivation can be found in the 
Appendix of the third of these references. The kinetic energy sum can be done analytically for calculations performed 
in a Jacobi basis.) This reordered BH equation has the form 



H^ff = H + HQ ^ ^^ OH ^ 



(2) 



E-QT 



E-QH^ E-TQ 
where the bare H is the sum of the relative kinetic energy and a two-body interaction 

iJ = - y (T,,, -I- y„) , with T„ = ilViP^L. (3) 

This effective interaction is to be evaluated between a finite basis of Slater determinants |a) G P, which is equivalent 
to evaluating the Hamiltonian 

H-lf =T- T^T + V + V ^ -QH ^^ ^"^^ 

between the states 

|S>-^H (5) 

By summing QT to all orders, the proper behavior at large r can be built in, which then allows b to be adjusted, 
without affecting the long-wavelength properties of the wave function. Fig. [l] from Ref. [B], shows that the resulting 



4 




I \ I \ I I I \ 

0.4 0.6 0.8 1.0 1.2 

b(fm) 



FIG. 1: (Color online) Deuteron ground-state convergence for "bare" calculations in small P-spaces, which omit all effects due 
to the multiple scattering of V in Q. The three curves on the upper right were calculated from the standard BH equation, 
which identifies the bare interaction as P{T + V)P. These calculations fail to bind the deuteron, even with Ap = 10, for all 
values of the HO size parameter b: the P-space estimate for V is poor if b is much above Ifm, while the estimate for T is poor 
if b is below that value. The lower four curves were evaluated for the bare interaction of the reordered BH given by Eq. ([2| , 
which incorporates the long-range effects of QT to all orders, building in the correct asymptotic form of the wave function. This 
allows one to reduce b to small values, pulling most of the effects of V into P, without distorting the long-distance behavior 
of the wave function or, therefore, the estimate for T. Rather remarkably, this bare calculation reproduces the correct binding 
energy for P spaces as small as Ap—6. That is, by the combination of the summation of QT to all orders and the adjustment 
of b to an optimal value characteristic of the hard core radius of vis, the effective interaction contribution can be driven to such 
small values that it can be ignored. 



decoupling of the long- and short- wavelength physics can greatly improve convergence: a "bare" 6 hcu calculation that 
neglects all contributions of QV gives an excellent binding energy. This decoupling of QV and QT is also important 
in finding a systematic expansion for TT^^-f . 

This reorganization produces an H'^-^^ with three terms operating between HO Slater determinants, 

(«|T— ^|/3) = (a|— ^T|/3) ""^^'^ {a\T\(3) 
^ ^ E ~QT^' ' ^ ^E-TQ \ I IH/ 

The ladder properties of QT make E/(E — QT) the identity operator except when it acts on an \a) with energy Kpfiuj 
or (Ap — \)fiuj. These are called the edge states. For nonedge states, the new grouping of terms in H'^^^ reduces to 
the expressions on the right-hand side of Eq. [g] the conventional components of TT^^^ . Thus the summation over QT 
alters only a subset of the matrix elements of H'^^^ , while leaving other states unaffected. 

Figure |2] shows the extended tail of the relative two-particle wave function that is induced by E/{E — QT) acting 
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FIG. 2: (Color online) A comparison of the radial wave functions for the HO state \nl) (dashed) and for the extended state 
{E/E — QT)\nl) (solid), for (n, I) = (6, 0) in a Ap = 10 deuteron calculation. The extended tail of the latter is apparent. Note 
that the normalization of the extended state has been adjusted to match that of \nl) at r=Q, in order to show that the shapes 
differ only at large r. Thus a depletion of the extended state at small r is not apparent in this figure. 



on an edge HO state [6]. As will become apparent from later expressions, this tail has the proper exponential fall-off. 



(7) 



where k = ^2\E\/hu) and r = |ri — r2|/\/2& is the dimensionless Jacobi coordinate, not the Gaussian tail of the HO. 
At small r the wave function is basically unchanged (apart from normalization). 



III. THE HOBET EFFECTIVE INTERACTION 



Contact-gradient expansions are used in approaches like EFT to correct for the exclusion of short-range (high- 
momentum) interactions. The most general scalar interaction is constructed, consistent with Hermiticity, parity 
conservation, and time-reversal invariance, as an expansion in the momentum. Such an interaction for the two- 
nucleon system, expanded to order N'^LO (or up to six gradients), is shown in Table |l] (Later these operators will be 
shghtly modified for HOBET.) 

The "data" for testing such an expansion for HOBET are deuteron matrix elements {a\P{H'^^^ — H)P\(i) evaluated 
as in Refs. [SJ|B] for wig. I take an ^fiuj P-space (Ap — 8). The evolution of the matrix elements will be followed as 
contributions from scattering in Q are integrated out progressively, starting with the highest energy contributions. To 
accomplish this, the contribution to W^^^ coming from excitations in Q up to a scale A > Ap is defined as H'^-f^{A), 
obtained by explicitly summing over all states in Q up to that scale: 

1 A 

H'ff{A) = H + H QaH Qa= ^ \a){a\ Qa, = 0. (8) 

Thus H^ff = H^f^A oo) and H^^^iAp) = H. The quantity 

A(A) . H^ff - H^ffiA) . H^QH - H^^Q.H (9) 



represents the contributions to H'^ff involving excitations in Q above the scale A. For A >> Ap, one expects A (A) 
to be small and well represented by a LO interaction. As A runs to values closer to Ap, one would expect to find that 
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TABLE I: Contact-gradient expansion for relative-coordinate two-particle matrix elements. Here Dm— (V ® V)2Af, D1= 

[{G{l)®a(2))2®D\o, Fli= (V ® D^)3M, Fl,= [(o-(l) ® cr(2))2 ® i^^'liM, Gt,= (D^ ® D^)iM, Gli= [{o{l) ® (j{2))2® G^m , 
and the scalar product of tensor operators is defined as A'' ■ B'' = X^jJ^Zl j(— 1)*^^M^-Af • 



Transitions 


LO 


NLO 


NNLO 


N^LO 


or ^5o ^ """So 




^'ioCV' 5(r)-f5(r) V2) 










a^^o(5(r)D0 + D0 5{r)) 




4?Lo(^' '5(r) DO + DO S{r) V^)^ 


or ^Dj ^ ^Dj 




















or ^ 3p^ 






<i^Lo(W -^W V + V -Sir) V^V) 


"^"lo ^ 








<^Lo(^ -5^ Fi + fi .5(r) V) 




or ^Fj ^ ^Fj 











NLO, NNLO, N^^LO, .... contributions become successively more important. If one could formulate some expansion 
that continues to accurately reproduce the various matrix elements of A (A) as A — > Ap, then a successful expansion 
for the HOBET effective interaction A(Ap) = H'^^^ — H would be in hand. 

Figure [S^ is a plot of A(A) for the 15 '^S'l matrix elements in the chosen P-space. For typical matrix elements 
A(Ap) = H^'-ff - H ^ -12 MeV - a great deal of the deuteron binding comes from the Q-space. Five of the matrix 
elements involve bra or ket edge states. The evolution of these contributions with A appears to be less regular than 
is observed for nonedge-state matrix elements. 

One can test whether the results shown in Fig. can be reproduced in a contact-gradient expansion. At each A 
the coefficients a|;U(A), a'^2o(^)' ^^c., would be determined from the lowest-energy "data," those matrix elements 
(q;|A(A)|/3) carrying the fewest HO quanta. Thus, in LO, o^j^q^K) would be determined from the (n',n) = (1,1) 
matrix element. The remaining 14 P-space matrix elements are then predicted, not fit; in NNLO four coefficients 
would be determined from the (1,1), (1,2), (1,3), and (2,2) matrix elements, and eleven predicted. Figures [sja-d 
show the residuals - the differences between the predicted and calculated matrix elements. For successive LO, NLO, 
and NNLO calculations, the scale at which residuals in A are significant, say greater than 10 keV, is brought down 
successively, e.g., from an initial ^ lOOTjo-i, to ^ BOTiw (LO), to ~ 30?ia; (NLO), and finally to ~ 20?iw (NNLO), 
except for matrix elements involving edge states. There the improvement is not significant, with noticeable deviations 
remaining at ~ 100?iw even at NNLO. This irregularity indicates a flaw in the underlying physics of this approach 
- specifically the use of a short-range expansion for H'^^^ when important contributions to H'^^f are coming from 
long-range interactions in Q. So this must be fixed. 



A. The contact-gradient expansion for HOBET 

The gradient with respect to the dimensionless coordinate r = {ri — r-2)lh\f2 is denoted by ^ . The coefficients 
a-LOi O'NLO, ■■■ in Table [l] then carry the dimensions of MeV. 

The contact-gradient expansion defined in Table |l] is that commonly used in plane- wave bases, where one expands 
around k = with 



^ exp ik ■ r 



, =0. (10) 

fe=0 



7 






FIG. 3: In a) the contributions to H"" — H from excitations in Q above A are plotted for a calculation with Ap = 8 and 
6=1.7 fm. Each line describes the running of one of the 15 independent P-space matrix elements {n'l' = OlH"^-^ — H\nl = 0), 
1 < n < n' < 5. Ten of the matrix elements are between nonedge states (solid), four connect the n' = 5 edge state to the 
n=l,2,3,4 nonedge states (dashed), and one is the diagonal n' = n — 5 cdgo-odgc case (dot dashed). b)-c) show the residuals 
for naive LO, NLO, NNLO, and N^LO fits (see text), f) shows the RMS deviation for the set of P-space matrix elements. The 
expected systematic improvement with increasing order is apparent only for matrix elements between nonedge states. 
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HOBET begins with a lowest-energy Is Gaussian wave packet with a characteristic momentum 1/b. An analogous 
definition of gradients such that 

^Vi.W = (11) 
is obtained by redefining each operator appearing in Table |l] by 

0^0 = e'-'/20e'-'/2. (12) 

The gradients appearing in the operators of Table |l] then act on polynomials in r. This leads to two attractive 
properties. First is the removal of operator mixing. Once is fixed in LO to the (n', n) = (1, 1) matrix element, 
this quantity remains fixed in NLO, NNLO, etc. Higher-order terms make no contributions to this matrix element. 
Similarly, aMLO, once fixed to the (1,2) matrix element, is unchanged in NNLO. That is, the NLO results contain the 
LO results, and so on. Second, this definition gives the HOBET effective interaction a simple dependence on nodal 
quantum numbers, 

"^^--4(71-1) - 16(n- l)(n- 2). (13) 

(The Appendix describes this expansion in some detail.) In each channel, this dependence agrees with the plane- 
wave result in lowest contributing order, but otherwise differs in terms of relative order 1/n. This HO form of the 
contact-gradient expansion is connected with standard Talmi integrals [:8l, generalized for nonlocal potentials, e.g., 

(14) 

Jo Jo Jo Jo 

and so on. 



aLo - / e ''^ [V{ri,r2)]e ''^r(r^dridr2 

Jo Jo 

/•OO /"OO /"OO /"OO 

aNLO-- / / e"''! [r?y(ri,r2)] e-'^^V^r^dridra = / / e~'^^V{ri,r2)rl] e''^^^rlrldridr2 
Jo Jo Jo Jo 



B. Identifying terms with the contact-gradient expansion 

The next question is the association of the operators in Table [l] with an appropriate set of terms in H'^^^ — H, so 
that the difficulties apparent in Fig. [sjare avoided. The reorganized BH equation of Eq. ([2| 



E-TQ 
E 

E-TQ 



t-tI^t + v + v^^qv 

T - T—T + V+ V O' 
E ^ 

i=LO,NLO.... 



E 



E-QT 
E 

E-QT 



(15) 



isolates V{E — QH)~^QV , a term that is sandwiched between short-range operators that scatter to high-energy 
states: one anticipates this term can be successfully represented by a short-range expansion like the contact-gradient 
expansion. This identification is made here and tested later in this paper. 

This reorganization only affects the edge-state matrix elements, clearly. As the process of fitting coefficients uses 
matrix elements of low (n', n), none of which involves edge states, the coefficients are unchanged. But every matrix 
element involving edge states now includes the effects of rescattering by QT to all orders. Thus a procedure for 
evaluating these matrix elements is needed. 

C. Matrix element evaluation 



There are several alternatives for evaluating Eq. (151 for edge states. One of these exploits the tri-diagonal form of 



QT for the deuteron. If \nl) is an edge state in P then 

E^Q^ l^ = 1^0+ E-QT ^^^"" = \nl) + ^nin + I + 1/2) _^^[^_^ \n + 1 /) (16) 
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'2\E\ 

where E < Q for a bound state. The dimensionless parameter n — \ -j-j- depends on the ratio of the binding energy 
the HO ener, 

-QT\n + l I) 



\E\ to the HO energy scale. Note that the second vector on the right in Eq. (16 1 hes entirely in Q. Now 
2 



2 

flLO 

2 

TlLO 

2 

TlLO 



QT\n + 2 I) 
QT\n + i I) 
QT\n + Al) 



(2n + Z + 3/2)|?i + 1 I) + y/{n + l){n + l + i/2)\n + 2 I) 
^{n+\){n + l + 'i/2)\n + U) + (2n + / + 7/2)|n + 2l) + ^ {n + 2){n + I + b/2)\n + 30 
y/{n + 2){n + l + b/2)\n + 2 1) + {2n + l+ ll/2)|n + 3 + {n + ■i){n + I + 7 /2)\n + A I) 



(17) 



So the operator 2QT jYiuj in the basis {|n + i Z) , * = 1, 2, ... } has the form 

/ ai /3i \ 

/?i a2 /32 



fii as /^a 
/33 a4 



(18) 



where 



cti = ai{n,l) =^2n + 2i + l- 1/2, /3j = /3i(n, ?) = ^/lJ^Tij{(rlTiTTTT/fj . 



(19) 



As is well known, if this representation of the operator 2QT/'huj is truncated after k steps, the 2fc — 1 nonzero 
coefficients {a^. Pi} determine the 2fc — 1 operator moments of the starting vector |n + 1 I), 



{n+ll\ 



(2QTV 
\ huj J 



11), i = l,....,2fc- 1 



(20) 



A standard formula exists ^ for the moments expansion of the Green's function acting on the first vector |n + 1 I) 
of such a tri-diagonal matrix, allowing us to write 



^yn{n + I + 1/2)- 



1 



K-2 _ 2QT 



|n+l I) = gi{-K'^;n,l)\n+l I) + g2{-K'^;n,l)\n + 2 I) + g3{-K'^;n,l)\n + 3 I) •{ (21) 



The coefficients {gi} can be obtained from an auxiliary set of continued fractions {g[} that are determined by downward 
recursion 



g'^.{~K^\n,l) = 



1 



-k2 - ak{n,l) 



1 



9i-i{ K^',n,l) — ( i\ Q ( nn I ( 1 i\' * — •••■2 

From these continued fractions the needed coefficients can be computed from the algebraic relations 

gi{-K^;n,l) = y'nin + I + l/2)g[{-K'^;n,l) 
gi{-K^;n,l) = g,,_i(-K^; n, /)/3j_i5-(-k^; n, ?), i = 2, 



(22) 



(23) 



Defining gQ{~K^\n,l) = 1 it follows 



E 



E-QT' 



k — )-oo 

l> = gi{—K^;n,l)\n + i I), edge state 

1=0 

— \n I otherwise 



(24) 



where it is understood that k is made large enough so that the moments expansion for the Green's function is accurate 
throughout the region in coordinate space where E/{E ~ QT)\n I) is needed. Note that the ffist line of Eq. (24) can 
be viewed as the general result if one defines 



gi{—K^; n,l) = 0, i — 1, ...k, if \n I) is not an edge state. 



(25) 
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(For A > 3 one would be treating the 3(^-l)-diniensional HO, with the role of the spherical harmonics replaced by 
the corresponding hyperspherical harmonics.) Eq. (24) can now be used to evaluate the various terms in Eq. (15 1. 



Matrix elements for the contact-gradient operators: The matrix elements have the general form 



{n'l'l 



E 



:0- 



E 



E-TQ E-QT' 



^0 = X! 9]{-i^'^;n' ,l')gi{-K^ 



n,l){n' +j l\0\n + i I) 



(26) 



where O is formed from gradients acting on the bra and ket, evaluated at r—0. The general matrix element (any 
partial wave) is worked out in the Appendix. For example, one needs for S'-wave channels the relation 



(n-iy. 1 
{n — 1 — p)\ n 



r(n + 1/2) 1^/' 
(n-1)! 



(27) 



from which it follows 
{n'{l' ^0 5*= 1)J= 1| 



E 



E-TQ 



E 



E 



E-QT 



V gj(-K2; n', I' = 0)U-^^; n, I = 0) 

id=o 



Q3Sl,i 

i=LO,...,N^LO 

Tn' + j + l/2Tn + i+l/2 



|n(Z = 5*= 1)J= 1) 



,3S1 
''LO 



{n' + i - l)\{n + i - 1)! 

+ 16 [{n' + J - l){n + I - l)a^^l^ll + + J - +]-2) + {n + i- l){n + i - 2)) a^f^'^^} 
-Qa[{7i'+] - l){n + I - l)((n' + J -2) + {n + i- 2))a%\^^ 



4((n' + j-l) + (n + i-l))a3f2o 



+ ((n' + j - l)(n' + i- 2){n' + j - 3) + (n + i ~ l)(n + i - 2)(n + i - 3)) 



351,60 I 



(28) 



In the case of nonedge states, gi = Q except for the case of go = 1- Thus it is apparent that the net consequence of 
the rearrangement of the BH equation and the identification of the contact-gradient expansion with V{E — QH)~^QV, 
is effectively a renormalization of the coefficients of that expansion for the edge HO states. That renormalization is 
governed by = 2\E\/hLLj, e.g., 



aLo{n\ I', n, I) a'^^oi^' ^' ^ = a,Lo{n', I', n, I) ^ 9j{-n^\ n' , l')gi{-n'^; n, I) 



r(n'+j + l/2)r(n + i + l/2) 
r(n' + l/2)r(n+l/2) 



1/2 



(n'-l)!(n-l)! 



(n' +j - iy.{n + i - 1)! 



1/2 



(29) 



This renormalization is large, typically a reduction in strength by a factor of 2-4, for \E\—2.22A MeV, and also 
remains substantial for more deeply bound systems, as will be illustrated later. (The binding energy for this purpose 
is defined relative to the lowest particle breakup channel, the first extended state.) The effects encoded into \a) by 
summing QT to all orders are nontrivial: they depend on a nonperturbative strong interaction parameter \E\ as well 
as QT, and they alter effective matrix elements of the strong potential. For a given choice of Ap, the renormalization 
depends on a single parameter, 2\E\/?iuj, not on \E\ or b separately. In the plane-wave limit b oo, this parameter 
is driven to oo, so that a'j^Q ulo- No renormalization is required in this limit. The dependence on \E\ is discussed 
in more detail later, including its connection to the state-dependence inherent in effective theory. 

Matrix elements of the relative kinetic energy: The relative kinetic energy operator couples P and Q via strong matrix 
elements that grow as n. As Ref. [6., discusses, this coupling causes difficulties with perturbative expansions in H 
even in the case of P spaces that contain almost all of the wave function (e.g., Ap ^ 70). There is always a portion 
of the wave function tail at large r that is nonperturbative, involving matrix elements of T that exceed KpTiLo/2. 
The kinetic energy contribution is 



{a\T + T 



1 



E-QT 



QT\P) = {a\T-T^T\l3) = {a\T\P) 



{a\Tl 



(30) 



where the last two terms show that the transformation to states \a) — E/{E — QT)\a) reduces the calculation of the 
rescattering to that of a bare matrix element. It follows from this expression 



in' l\T + T—^—QT\n I) = {n' l\T\n I) + ^5„'„x/"(" + ' + 1/2)5i(-k'; n, /). 
hi — Qi z 



(31) 
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Thus, rescattering via QT alters the diagonal matrix element of the effective interaction for edge states, as determined 
by 5i(-K^;n,0- 

Matrix elements of the hare potential: The P-space matrix element of V becomes (a|y|/3) which, as is illustrated in 
Fig. [2] involves an integral over a wave function that, apart from normalization, differs from the HO only in the tail, 
where the potential is weak. It can be evaluated by generating the wave functions \a) and as HO expansions. 



V 



^ |n + i l)g,{-K^;n,l) 



.4=0 



gj{-K^-n' ,l')g,{-K^;n,l){n' + j l'\ V \n + il). 

(32) 

though the alternative Green's function expression, discussed below, is simpler. 



Use of the free Green's function: An alternative to an expansion in an HO basis is generation of \a) with the free 
(modified Helmholtz) Green's function. For any P-space state \n I), 

{E ~ QT)\a) = E\a) ^ {E ~ T)\a) ^ E\a) ~ PT\a) (33) 

That is, both E — QT and E~T project \a) back into the P-space. The free Green's function equation can be written 

E 



{E-T)\a) 



P 



E-T- 



P- 



1 

E~T' 



E-QT 
-I -1 



P\a) 



-P 



(34) 



Either of the driving terms on the right-hand side is easy to manipulate. The second expression requires inversion of 
a P-space matrix, one most easily calculated in momentum space, as the HO is its own Fourier transform and as the 
resulting momentum-space integrals can be done in closed form. This form was used in the three-body calculations 
ofRef. [TD]. 

Here I will use the first expression a^bove, rewriting the right-hand-side driving term in terms of Iq^^^^t^^), 



p 



E-T- 



E 



E-QT 



P\a) - 



flOJ 



-n'' -{2n + l-\/2)-M' 



-K^; n, l)^n{n + / + 1/2)^ \n I mi > 



V("- - l)in + l - l/2)|n - 1 I mi) - ^/nin + I + l/2)P\n + 1 / mi) 



(35) 



where the driving term has been kept general, valid for either edge or nonedge states: the latter can be a helpful 
numerical check, verifying that a HO wave function is obtained, for such cases, from the expression below. For an 
edge state, gi is nonzero and P\n + 1 I) = 0; for a nonedge state, gi = and P=l. Labeling the corresponding edge 
state as |a„im,). 



{r\anim,) 



(Pr 



1 



47r|r — f 

^Im (^r) 



1 f 

-r Ii+i/2{^r) / d^r'{r')'^^Ki+^/2{Kr') (r"|a„,„, > 

V Jr 



+ ^ AVi/2(«r-) / dV-"(/)3/'/z+i/2(Kr') (f'|a„,„, 



> 



(36) 



where / and K denote the standard modified Bessel functions. By expressing the HO radial wave functions in terms 
of the underlying Laguerre polynomials and integrating the polynomials term by term, alternative expressions are 
obtained for the various quantities previously expressed as expansions in the gi. This is detailed in the Appendix. 
One finds, for example. 



0|5, 



nlmi I 



(n-l)!r(n-H/2) ^ 



27r 



fe=0 



fc!(n-fc)!r(fc-^3/2) 



{n - k){K^ + 371 - 3/2 - fc + gi(-K^; n, 0)V"(" + 1/2)) + P{n + 1, / = 0]n(n + 1/2) 
-^2 K r(fc + 3/2) iPi[fc-H3/2;3/2;KV2] + iPi [fc + 1; 1/2; kV2] 



(37) 
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where iFi is Kummer's confluent hypergeometric function, and where P[n + 1, Z = 0] = 1 if |n + 1 /) is in P, and 
otherwise. Similar expressions can be derived to handle all of the operators O appearing in the contact-gradient 
expansion (see the Appendix). 



D. Numerical Tests 



In this subsection channel-by-channel N'^LO results are presented for H'^^^ based on Eqs. (12 1 and (15 1, which 
isolate a short-range operator that plausibly can be accurately and systematically expanded via contact-gradient 
operators. 

For Ap — 8, the fitting procedure determines all N"^LO coefficients from nonedge matrix elements, leaving all edge 
matrix elements and a substantial set of nonedge matrix elements unconstrained. Thus one can use these matrix 
elements to test whether the expansion systematically accounts for the "data," the set of numerically generated vis 
matrix elements of H'^f^. One test is the running of the results as a function of A: a systematic progression through 
LO, NLO, etc., operators should be observed as A is lowered to the SM scale. A second test is the "Lepage plot" 
|llj . which displays residual errors in matrix elements: if the improvement is systematic, these residual errors should 
reflect the nodal-quantum-number dependence of the operators that would correct these results, in next order. 



Eq. ( 15 1 includes "bare" terms ~ the matrix elements {a\T\l3) and (S|F|/3) - and a term involving repeated scattering 
by H in Q, but sandwiched between the short-range operator QV. To test the dependence on A, the rescattering 
term is decomposed in the manner ofEq. ([9|, 



E 

Aqt(A) = 



E-TQ 



E 



E-QT' 



to isolate the contribution of scattering above the scale A. Aqt(A) is evaluated numerically for vis at each required 
A. The long-wavelength summation is always done to all orders - the running with A thus reflects the behavior of the 
short-range piece, V{E — QH)^^QV . The full P-space effective interaction is obtained as A ^ Ap. 

As outlined before, coefficients are fitting to the longest wavelength information. For example, in S channels, a^o 
is fixed to the {n',n) = (1,1) matrix element; the absence of operator mixing then guarantees this coefficient remains 
fixed, as higher order terms are evaluated. The single ajvLO coefficient is fixed to (2,1) (or equivalently (1,2)) ; a^jvLO 
and ajvAfio determined from (2,2) and (3,1); and finally o^^-^^q and a^si^p ^'-^ fixed to (3,2) and (4,1). So at 
N'^LO there are a total of 6 parameters. This procedure is repeated for a series of A ranging from 140 to Ap—8. 
The results in each order, and the improvement order by order, are thus obtained as a function of A. P contains 15 
independent matrix elements in the ^Si — '^Si channel, nine of which play no role in the fitting: these test whether 
the improvement is systematic. 

Figures H and [s] show the results for ^Si —^Si and ^5*0 — ^5*0. Panel a) shows the evolution of the matrix elements 
{a\AQT{Aj\f3) for each of the 15 independent matrix elements. Matrix elements involving only nonedge states, a single 
edge state, or two edge states are denoted by solid, dashed, and dash-dotted lines, respectively. Progressively more 
binding is recovered as A — > Ap. In the — '^S'l case, the contribution at Ap is ~ 12-14 MeV for nonedge matrix 
elements, ~ 7-8 MeV for matrix elements with one edge state, and ^ 5 MeV for the {n = 51 = 0|AQT(Ap)|n = 5^ = 0) 
double-edge matrix element. 

Panels b)-e) show the residuals - the difference between the matrix elements of Aqt(Ap) and those of the contact- 



gradient potential of Eq. (15) - from LO through N'^LO. The trajectories correspond to the unconstrained matrix 
elements (14 in LO, 9 in N^LO): the fitted matrix elements produce the horizontal line at 0. Unlike the naive 
approach in Fig. |3j the improvement is now systematic in all matrix elements. In the ^Si — "^Si channel, a LO 
treatment effectively removes all contributions in Q above A ~ 60; NLO lowers this scale to ~ 40, and NNLO is ^ 
20. The magnitude of N'^LO residuals at Ap is typically < 2 keV - the entire effective interaction can be represented 



by Eq. (151 to an accuracy of about 0.01%. Panel f) shows the root mean square (RMS) deviation among the 
unconstrained matrix elements, and the rapid order-by-ordcr improvement. 

The pattern repeats in the ^5*0 — channel, where the convergence (in terms of the size of the residuals) is 
somewhat faster. The N'^LO RMS deviation among the unconstrained matrix elements at Ap is ^ 0.5 keV. 

The remaining positive-parity channels that can be constrained at N^LO are given in Figs. [6| through |ll[ '^5'i — ^Di 
(leading order contribution NLO); ^D2-^D2, ^Di-^Di, ^D2-^D2, and -^Dg (NNLO); and^Dg-^g (N^LO). 

Table H gives the resulting fitted couplings at N^LO for all contributing channels, along with numerical results for 
the root-mean-square Q-space contributions to Aqt(Ap) and the root- mean-square residuals (the deviation between 
the contact-gradient prediction and Aqt(Ap) for the remaining unconstrained effective interactions matrix elements). 
The quality of the agreement found in the ^Sq and 5*1 channels is generally typical ~ residuals at the few kilovolt 
level - though there are some exceptions and some general patterns that emerge. One of these is the tendency of the 
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> 



C) 





FIG. 4: As in Fig. [s] but for the QT-summed reordering of W^^^ . The contributions to the effective interaction from excitations 
in Q above A, denoted Aqt(A) in the text, are plotted. Each line gives the running of a P-space matrix element. b)-e) show the 
residuals for LO, NLO, NNLO, and N'^LO fits (see text), f) shows the RMS deviation for the set of P-space matrix elements. 
The improvement with increasing order is systematic and rapid: at N^LO the RMS deviation for unconstrained matrix elements 
as A ^ Ap is about 3 keV. That is, the entire effective interaction is reproduced to a few parts in 10*. 
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FIG. 5: As in Fig. [4] but for the ^So channel. The N'^LO results are seen to reproduce the entire effective interaction to the 
accuracy of about a keV, or one part in 10^. 
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FIG. 6: As in Figs. |4]and|5l but for the ^Si - ^Di channel. 
As with the cases described before, the N'^LO results remain 
accurate at the few keV level, as the integration is brought down 
to the shell-model scale, A ^ Ap. 
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FIG. 7: As in Fig. [i] but for the ^Di channel. 



triplet channels with spin and angular momentum aligned (^S\ — ^S\, — ^P2, ^Ds — "^Ds, and -^^4 — ^FA) to exhibit 
larger residuals than the remaining S,P,D and F channels, respectively. The '^£'3 — "^D^, which has contributions at 
NNLO and N'^LO, stands out as the most difficult channel, with a residual of 122 keV, one to two orders of magnitude 
greater than the typical scale of N'^LO residuals. 

Figures [12] through [18] show the convergence for the various channels involving odd-parity states and contributing 
through N^O: ^Pi -^1, ^Pj - ^Pj, ^P2 - ^^2, ^-^3 - ^^3, and ^Fj - ^Fj. While the spin-ahgned channels show 
slightly large residuals, overall the RMS errors at N'^LO are at the one-to-few keV level. Thus a simple and essentially 
exact representation for the effective interaction exists. 



Expansion parameters, naturalness: The approach followed here differs from EFT, where the formalism is based on 
an explicit expansion parameter, the ratio of the momentum to a momentum cutoff. The input into the present 
calculation is a set of numerical matrix elements of an iterated, nonrelativistic potential operating in Q. Potentials 
like 1118 are also effectively regulated at small r by some assumed form, e.g., a Gaussian, matched smoothly to the 
region in r that is constrained by scattering data. Thus there are no singular potentials iterating in Q. 

Intuitively it is clear that the convergence apparent in Table jll] is connected with the range of hard-core interactions 
(once edge states are transformed by summing T). A handwaving argument can be made by assuming rescattering 
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in Q effectively generates a potential of the form 



2 / 2 



where ri2 = l^i — r2|. This ansatz is local, so there is some arbitrariness in mapping it onto contact-gradient expansion 
coefficients, which correspond to the most general nonlocal potential. But a sensible prescription is to equate terms 
with equivalent powers of r^, in the bra and ket, when taking HO matrix of this potential. Then one finds, for S-wave 
channels 



a(m', m) = d 



S,2m'2m 



1 



(2m' + 2TO+l)!! 



4™'+™TO'!m! (2m' + 1)!!(2to + 1)!! 



^^0 





3/2 






r a" 1 


a2 + 262 




a2 + 262 



m +m 



(38) 



where h is the oscillator parameter. (The notation is such that, e.g., a{m' = 3,m = 0) = 



S,60 



) The last term is 



thus the expansion parameter: if the range of the hard-core physics residing in Q is small compared to the natural 
nuclear size scale 6, then each additional order in the expansion should be suppressed by ~ (a/5) 2. 

One can use this crude ansatz to assess whether the convergence shown in Table is natural, or within expectations. 
The LO and NLO ^5*0 — ^S'o results effectively determine Vb and a; thus the strengths of four NNLO and N'^LO 
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FIG. 9: As in Fig. [t] and (s) but for the '^Do, channel. This 
"stretched" configuration generates much larger residuals than 
the other 1 — 2 channels. Consequently a calculation to N'^LO 
would be needed to reduce typical matrix element errors to ~ 
10 keV, in the limit A — > Ap. 
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FIG. 10: As in Fig. ^ but for the ^732 channel. 



potentials can be predicted relative to that of qlo and unlo- The predicted hierarchy ulo ■ o-nlo ■ O'j^NLO 

„40 . 42 . 60 „f 



1 : 6.3 X 10"^ : 6.7 x 10"^ : 2.0 x lO"'^ : 3.0 x 10"^ : 4.2 x 10"* 

matches the relative strengths of the couplings in the table quite well, 

1 : 6.3 X 10"^ : 6.4 x 10"^ : 3.9 x 10"^ : 2.1 x 10"^ : 5.7 x 10"^ 

including qualitatively reproducing the ratios of the two NNLO and two N'^LO coefficients. The parameters derived 
from a^o and a^LO are a ~ 0.39 fm and Vq ~ —1.5 GeV. In the ^Sq channel the bare Argonne ?;i8 potential at small 
r can be approximated by a Gaussian with a 0.33 fm and Vq ~ 3.0 GeV. So again the crude estimates of range 
and even the strength are not unreasonable. [Note that the signs of the two VqS are correct - the P-space lacks the 
appropriate short-range repulsion and thus samples the iterated bare potential at small r, a contribution that then 
must be subtracted off when H'^-^^ is evaluated.] 

A similar exercise in the '^Si channel yields the predicted hierarchy 



1 : 2.2 X IQ-^ : 8.3 x 10"'' : 2.5 x 10"^ : 13.1 x 10"*^ : 1.9 x 10" 
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FIG. 11: As in Fig. [i) but for the ^D^ - ^Ga channel. The N*LO contribution is also shown. 



which compares with the coupl: 



ing ratios calculated from Table [lT| 
1 : 2.2 X 10"^ : 8.8 x lO"'^ : 1.9 x 10"* : 9.3 x 10"^ : 2.3 x 10"^ 



The convergence is very regular but slower: in this case the effective Gaussian parameter needed to describe these 
trend is a ~ 0.75 fm. The overall strength, Vq ^ —0.42, differs substantially from that found for the ^Sq channel, 
though the underlying vis potentials for — ^Si and ^5*0 — "^^o scattering are quite similar (see Fig. 19). 



The '^Si — "^Si behavior is similar to that found in the other spin-aligned channels, such as ^D^ — ^D^ and — "^^2, 
where the scattering in Q includes contributions from the tensor force. The tensor force contributes to the LO s-wave 
coupling through intermediate Z?-states in Q, e.g., 



{n'l' = Q\VsDQ\n"l 



2\QVsD\nl^Q), 



as the product of two tensor operators has an s-wave piece. The radial dependence of Vsd for wig, shown in Fig. [19] 
is significantly more extended than in central-force ^5*1 — ^S'l and ^5*0 — ^5*0 cases. This has the consequences that 

(1) the mean excitation energy {E) for ^S'l — ^Di will be lower (enhancing the importance of the tensor force) and 

(2) the P-space ^ Si\H'^^ ^ \^ Si) matrix element will reflect the extended range. 
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TABLE II: The effective interaction for LO through N'LO, with Ap = 8 and 6=1.7 fm. 
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* The appropriate LO, NLO, and NNLO interactions are obtained by truncating the table at the desired order, 
t An N^LO calculation in the 'P>3 - 'P3 channel yields a^^''^''=-2.510E-4 MeV and a^/^''^ = -7.550E-5 MeV, and reduces 
{Resid.)flMS to 22.3 keV; and in the 'D3 - '6*3 channel yields a^?^o = -2.141E-5 MeV and = 1.180E-5 MeV and 

reduces (Resid.) hms to 3.26 keV. 

Once this point is appreciated - that the effective expansion parameter are naturally channel-dependent because of 
effects hke the tensor force - the results shown in Table [iT] are very pleasing: 

• In each channel the deduced couplings a^o, clnlo, clnnlo, ^m^lOt-- evolve in a very orderly, or natural, 
fashion: one can reliably predict the size of the next omitted term. The convergence appears related to an 
effective range characterizing scattering in Q. 

• The convergence varies from channel to channel, but this variation reflects underlying physics, such as role of 
the tensor force, governing the channel's range. One does not find, nor perhaps should one expect to find, some 
single parameter p/A to characterize convergence independent of channel. 

• The convergence is very satisfactory in all channels: the measure used in Table [ll| (Resid.) raiSi is an exceedingly 
conservative one, as discussed below. But even by by this standard, in only one channel (^D^ — '^D^) do the RMS 
residual discrepancies among unconstrained matrix elements exceed ~ 10 keV. Given the arguments above, it 
is perfectly sensible to work to order NNLO in rapidly-converging channels like ^S'o — ^5*0 and N^LO in slowly 
converging channels like — ^D^. As noted in the table, at N'*LO the residual in the ■^£'3 — ^D^ channel is 
reduced to 22 keV. 

Convergence and the "Lepage" plot: The procedure often followed in an effective theory is to use information about 
the low-lying excitations to parameterize an effective Hamiltonian, which is then used to predict properties of other 
states near the ground state. In contrast, the goal here has been to characterize the entire effective interaction to high 
accuracy. As described below, the residual errors in the procedure are typically dominated by matrix elements with 
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FIG. 16: As in Fig. ^ but for the - channel. 



the largest n and n' , corresponding to minor components in the deuteron ground state, for example. The difference 
in the deuteron binding energy using exact matrix elements of H'^f^ versus using the N'^LO expansion is quite small 
(~ 40 eV). 

Order- by-order improvement should be governed by nodal quantum numbers. For example, in LO in the '^Si 
channel the omitted NLO term would be 



^^(n' + n-2) 



rK + i/2]r[n + i/2] 

{n' - l)\{n - l)\ 



1/2 , On^Sl 

"'^'■«'=-^^[(4n'-l)(4n 



l)]'/'(n' + n-2) 



(39) 



Thus the fractional error associated with the omission of the NLO terms relative to LO should to be linear in the 
sum of the nodal quantum numbers, if the expansion is capturing the correct physics. That is, the expected absolute 
(e.g., in keV) error for (n', rt)=(5,5) would be about 16 times that for (1,2). 

In higher orders this distinction between large and small n grows. At LO-I-NLO, the expected fractional errors in 
matrix elements from omitted NNLO terms would be quadratic in n and n': the explicit functional dependence is no 
longer simple as there are two NNLO operators, and one would not know a priori the relevant quadratic combination 
of n' and n governing the error. At NNLO the fractional error would be a cubic polynomial in n and n' . 

While beyond LO the expected fractional errors have a dependence on both n and n', it is still helpful to display 
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FIG. 17: The lowest contributing order to the IF3 channel is N^LO. Aqt(A) and the N^LO residuals for the five unconstrained 
matrix elements are shown. 



results as a 2D "Lepage plot" using n + n' - proportional to the average (p^) of bra and ket - as the variable. Such a 
plot makes clear whether improved fits in an effective theory are systematic - that is, due to a correct description of 
the underlying physics, not just additional parameters. The use of a single parameter, n + n', of course maps multiple 
matrix elements onto the same x coordinate, when the ET indicates this is a bit too simple beyond LO. Nevertheless, 
the right panel in Fig. [T9]still shows rather nicely that the nuclear effective interactions problem is a very well behaved 
effective theory. In LO the residual errors do map onto the single parameter n + n' to very good accuracy, and the 
residual error is linear. The steepening of the convergence with order is consistent with the expected progression from 
linear to quadratic to cubic behavior in nodal quantum numbers. By NNLO, errors in unconstrained matrix elements 
for small n + n' are tiny, compared to those with high n + n' . That is, the expansion converges most rapidly for 
matrix elements between long-wavelength states, as it should. However, improvement is substantial and systematic 
everywhere, including at the largest n + n' . 



IV. PROPERTIES AND ENERGY DEPENDENCE OF THE EFFECTIVE INTERACTION 



The results of the previous section demonstrate the existence of a simple systematic operator expansion for the 
HOBET effective interaction. Its behavior order-by-order and in the Lepage plot indicates that the short-wavelength 
physics is being efficiently captured in the associated operator coefficients. 

The error measure used in the N^LO fit is dominated by the absolute errors in matrix elements involving the 
highest nodal quantum numbers: these matrix elements are large even though they may not play a major role in 
determining low-lying eigenvalues. (It might have been better to use the fractional error in matrix elements, a measure 
that would be roughly independent of n' and n.) Other possible measures of error are the ground state energy; the 
first energy-moment of the effective interaction matrix (analogous to the mean eigenvalue in the SM); the fluctuation 
between neighboring eigenvalues of that matrix (analogous to the level spacing in the SM); and the overlap of the 
eigenfunctions of that matrix with the exact eigenfunctions (analogous to wave function overlaps in the SM). The 
N'^LO interaction in the coupled ^Si~^Di channel produces a ground state energy accurate to ^ 40 eV; a spectral first 
moment accurate to 1.81 keV; an RMS average deviation in the level spacing of 3.52 keV; and wave function overlaps 
that are unity to better than four significant digits. As rescattering in Q contributes ~ -10 MeV to eigenvalues, the 
accuracy of the N'^LO representation of the effective interaction is, by these spectral measures, on the order of 0.01%. 
As the best excited-state techniques in nuclear physics currently yield error bars of about 100 keV for the lightest 
nontrivial nuclei, this representation of the two-body effective interaction is effectively exact [II [T2]. 

The approach requires one to sum QT to all orders, producing a result that depends explicitly on \E\ - which 
in this context should be measured relative to the first breakup channel. While the associated effects increase with 
decreasing \E\, it will be shown later that the renormalization is substantial even for well-bound nuclear states. The 
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FIG. 19: (Color online) The left panel shows the radial dependence of the Argonne vis potential in the ^So — ^So, '^Si — ^Si, 
and ^Si — ^Di (tensor) channels. The last is clearly more extended. The right panel is a "Lepage plot" displaying fractional 
errors as a function of the order of the calculation, on log scales. The steepening of the slope with order is the sign of a well 
behaved, converging effective theory. 



deuteron is definitely not an extreme case. The effects are also sensitive to the choice of P, through b, which controls 
the mean momentum within P - a small b reduces the missing hard-core physics, but exacerbates the problems at 
long wavelengths, and conversely. Figure [T] suggests factor-of-two changes in the Q-space contribution to the deuteron 
binding energy can result from ~ 20% changes in b. At the outset, the dependence on and b seems like a difficulty 
for nuclear physics, as modest changes in these parameters alter predictions. 

One of the marvelous properties of the HO is that the QT sum can be done. The two effects discussed above turn 
out to be governed by a single parameter, k. The associated effects are nonperturbative in both QT and QV. In the 
case of QT an explicit sum to all orders is done. The effects are also implicitly nonperturbative in QV, because of 
the dependence on \E\. This is why the BH approach is so powerful: because \E\ is determined self-consistently, it is 
simple to incorporate this physics directly into the iterative process (which has been shown to converge very rapidly 
in the HOBET test cases A=2 and 3). When this is done, one finds that k affects results in three ways: 

• the rescattering of QT to all orders, T{E — QT)^^QT, is absorbed into a new "bare" matrix element {a\T\f3{K)); 

• the new "bare" matrix element {a{K)\V\P{K)) captures the effects of QT in all orders on the contribution 
first-order in V; and 

• the matrix elements of the short-range operators O, which contain all the multiple scattering of QV, are similarly 
modified, {a{K)\0\P{K)) . 

So far the discussion has focused on the problem of a single bound state of fixed binding energy \E\, the deuteron 
ground state. No discussion has occurred of expectations for problems in which multiple bound states, each with a 
different H'^ffdEl), might arise. But 1) the dependence of H'^ffdEl) on K arises already in the single-state case, 
which was not a priori obvious; and 2) state dependence (energy dependence in the case of BH) must arise in the case 
of multiple states, as this is the source of the required nonorthogonality of states when restricted to P, a requirement 
for a proper effective theory. So a question clearly arises about the connection between the explicit k dependence 
found for fixed \E\, and the additional energy dependence that might occur for a spectrum of states. 

Because other techniques, like Lee-Suzuki, have been used to address problem 2) , it is appropriate to first stress the 
relationship between k and the strong interaction parameters provided in Table ITTj The choice Ap=8 is helpful, as it 
shows there is no relation. Every short-range coefficient arising through order N^LO was determined from nonedge 
matrix elements: the fitting procedure matches the coefficients to the set of matrix elements with n' + n < 5, and 
there are no edge states satisfying this constraint. Nothing in the treatment of the strong interaction "knows" about 
edge states. This then makes clear how efficiently k captures the remaining missing physics. Without k one would 
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have, in the contact-gradient expansion to N^LO, a total of 78 poorly reproduced edge-state matrix elements, 10 of 
which would be S-state matrix elements with errors typically of several MeV. With k - a parameter nature (and the 
choice of b) determines - all of the 78 matrix elements are properly reproduced, consistent with the general ~ keV 
accuracy of the N'^LO description of H'^^-f . 

Suppose someone were to prefer an H'^f^ free of any dependence on again in the context of an isolated state 
of energy \E\. Could this be done? Yes, but at the cost of a cumbersome theory that obscures the remarkably simple 
physics behind the proper description of the edge state matrix elements. Suppose one wanted merely to fix the five 
'^Si—^Si edge state matrix elements, those where n' = 5 couples to n=l, 2, 3,4, and 5. One could introduce operators 
corresponding to the coefficients 

S,80 S,82 S,84 S,86 S,88 

^Ar*LO' ^N'^LO^ ^N'^LO' '^N''LO^ '^N^LO 

to correct these matrix elements. It is clear all five couplings would be needed - that's the price one would pay for 
mocking up long-range physics (a long series of high-order Talmi integrals) with a set of short-range operators of this 
sort. 

This would be a rather poorly motivated exercise: 

• The problems in these matrix elements have nothing to do with high-order generalized Talmi integrals of the 
strong potential, as was demonstrated in the previous section. 

• This approach does not "heal" the effective theory: the poor running of matrix elements would remain. There 
would be no systematic improvement, for all matrix elements, as a function of A, as one progresses from LO, to 
NLO, etc. The five parameters introduced above would remove the numerical discrepancies at Ap, but not fix 
the running as a function of A, even for just the edge-state matrix elements. 

• This approach amounts to parameter fitting, in contrast to the systematic improvement demonstrated in the 

S 80 

Lepage plot. The parameter a^4^Q introduced to fix the n = 1 to n = 5 matrix element will not properly 
correct the n — 2 to n — h matrix element, as the underlying physics has nothing to do with the rj-weighted 
Talmi integral of any potential. 

• If Ap is increased, the number of such edge-state matrix elements that will need to be corrected by the fictitious 
potential increases. This contrasts with the approach where \E\ is explicitly referenced: there the number of 
short-range coefficients needed to characterize Q will decrease (that is, the LO, NLO, ... expansion becomes more 
rapidly convergent) , while k remains the single parameter governing the renormalization of those coefficients for 
edge-state matrix elements. 

While these reasons are probably sufficient to discard any such notion of building a K-independent H'^f^ , consider 
now the consequences of changing b - which after all is an arbitrary choice. The short-range coefficients in Table [ll| 
will change: there is an underlying dependence on QV{fi2/b). This governs natural variations in the coefficients - 
one could estimate those variations based on some picture of the range of multiple scattering in Q, as was done in the 
"naturalness" discussion. But there would be additional changes in the ratios of edge to nonedge matrix elements, 
refiecting the changes in k. This would induce in any K-independent potential unnatural evolution in b. That is, the 
fake potential would look fake, as b is changed. 

The arguments above apply equally well to the case of the state-dependence associated with techniques like Lee- 
Suzuki. To an accuracy of about 95%, the K-dependence isolated in H'^^^ is also the state-dependence that one 
encounters when \E\is changed. This is a lovely result: the natural K-dependence that is already present in the case 
of a short-range expansion of H'^f^ for a fixed state, also gives us "for free" the BH state-dependence. The result 
is not at all surprising, physically: changes in \E\ will alter the balance between QT and QV , and that is precisely 
the physics that was disentangled by introducing k. Mathematically, it is also not surprising: changing b at fixed \E\ 
alters k, just as changes in |£'| for fixed b would. Thus all of the QT effects identified above, in considering a single 
state, must also arise when one considers spectral properties. 

This argument depends on showing that other, implicit energy dependence in H'^^^ is small compared in the explicit 
dependence captured in k. Such implicit dependence can reside in only one place, the fitted short-range coefficients. 



A. Energy Dependence 



The usual procedure for solving the BH equation. 
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involves steps to ensure self-consistency. As the energy appearing in the Green's function is the energy of the state 
being calculated, self-consistency requires iteration on this energy until convergence is achieved: an initial guess for E 
yields an H'^^^ [E) and thus an eigenvalue which then can be used in a new calculation of the interaction H'^f^{E'). 
This procedure is iterated until the eigenvalue coresponds to the energy used in calculating H'^^^. In practice, the 
convergence is achieved quite rapidly, typically after about five cycles. 

As the BH procedure produces a Hermitian H'^^^ , this energy dependence is essential in building into the formalism 
the correct relationship between the P-space and full-space wave functions, that the former are the restrictions of the 
latter (and thus cannot form an orthonormal set). This relationship allows the wave function to evolve smoothly to 
the exact result, in form and in normalization, as Ap — s- oo. 

Generally this energy dependence remains implicit because the BH equation is solved numerically: one obtains 
distinct sets of matrix elements {a\H'^^-f {Ei)\(3) for each state i, but the functional dependence on Ei is not immediately 
apparent. But that is not the case in the present treatment, where an analytic representation for the effective 
interaction has been obtained. 

While significant energy-dependent effects governed by k have been isolated, additional sources remain in the case 
of a spectrum of bound states. The identified energy-dependent terms are 

. {a\T+j^QT\P) = («|r|^(At)); 

• HE^V^^m = {a{K)\Vm^)); and 

. (a|^0^|/3) = {a{K)\Om^)) 
The implicit energy dependence not yet isolated resides in the coefficients of the contact-gradient expansion, 

• i^ivE^Qvm = HoiEm. 

To isolate this dependence, one must repeat the program that was executed for the deuteron ground state at a variety 
of energies, treating H^'ffdEl) as a function of \E\. The resulting variations in the extracted coefficients will then 
determine the size of the implicit energy dependence. Of course, all of the explicit energy dependence is treated as 
before, using the appropriate k. 

The simplest of the explicit terms is the "bare" kinetic energy 

{n'l\T\nT{K)) = {n'l\T+T—^—QT\nl) = {n'l\T\nl) + ^S^-nV^i^ + I + 1/2) n, I). 

-ti — LJl Z 

where effects only arise in the double-edge-state case. Two limits define the range of variation. As \E\ oo, gi 0, 
so the edge-state matrix element takes on its bare value, (2n + l — 1/2)Tilu/2. Similarly one can show gi{—K^; n,l) —> n 
as the binding energy \E\ approaches zero. Thus for small binding, the matrix element approaches {n + l — l/2)Tiuj/2. 
Thus the range is a broad one, nhuj/2, about 35 MeV for the parameters used in this paper. The behavior between 
these limits can be calculated. The results over 20 MeV in binding are shown in the upper left panel of Fig. [20] for 
S, P, and D states. One finds that even deeply bound (E=-20 MeV) states have very significant corrections due to 
QT: the scattering in Q reduces the edge-state kinetic energy matrix elements by (2-3) huj/2, which serves to lower 
the energy of the bound state. The kinetic energy decreases monitonically as \E\ — > 0. 

The second K-dependent term, the "bare" potential energy {a(n)\V\(3(K)) , is displayed over the same range in the 



lower left panel of Fig. 20 for the five '^S'l — ^Si edge-state matrix elements. These matrix elements are again quite 
sensitive to \E\, varying by 2-3 MeV over the 20 MeV range displayed in the figure. 

The third K-dependence is the renormalization of the contact-gradient coefficients for edge states, 

^''^^YQ^^B^W'" " M~^'-^^'^nM-^';nJ) {n'+j l'\ 0\n + tl) (40) 
Here O is fixed, while the explicit energy dependence carried by the gi (i.e., the effects of the interplay between QT 



and QV) is evaluated. The upper left panel in Fig. 20 gives the result for the diagonal edge-state matrix element, 
\n' I') — \n I) = |5 0). As has been seen in other cases, the reduction due to the QT — QV interplay is substantial 
throughout the illustrated 20 MeV range. Thus the large effects observed for the deuteron, a relatively weakly bound 
state, are in fact generic. But weakly bound states are more strongly affected, with the differences between the 
corrections for the double-edge states changing by a factor of nearly two between |i?|=20 MeV and \E\ ^ MeV. The 
results for single-edge-state matrix elements are similar, but the changes are smaller by a factor of two. 
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FIG. 20: Contributions to H^^^ with explicit energy depen- 
dence, for P defined by Ap = 8 and b — 1.7 fm. The upper left 
panel shows the diagonal "bare" kinetic energy term {a|r|a) 
for the edge states \a} = \n = 5 I = 0), \n = 4 I = 1), and 
|n = 4 Z = 2). The dots indicate the limiting values for very 
large and very small binding energies. The kinetic energy plot- 
ted is dimensionless, given in terms of fttj/2. The lower left panel 
gives the matrix elements of the bare potential V between ^5*1 
edge states, as a function of E. The upper right panel shows the 
evolution of the quantities a'i^Q{E\n' ,1' ,n,l) / aLoiji' ,1' ,n,l), 
a'f^i^Q{E;n' , I' ,n,l)/aNLo{n' ,1' ,n, I), etc., through N'^LO for 
the diagonal matrix element with \n — 5 I = 0). The general 
softening of such matrix elements is apparent, for small binding 
energy - repeated scattering by T through high-energy oscilla- 
tor states in Q spreads the wave function and thus reduces the 
effects of the strong potential at short range. This effect is car- 
ried by the edge states, because their renormalization is affected 
by the missing long-range physics. See the text for further dis- 
cussion. 



In doing these calculations, some care is needed in going to the limit of very small binding energies. One can show 
for edge states 



— / 2 , \ small K / \ 7 

g,i-K^;n,l) — > (-1)' 



T{n + l + l/2){n- 1 + 
T{n + I + 1/2 + i){n- 1)! 



1/2 



(41) 



If one uses this in Eq. (40 1 with n = 0, one finds that 

(r = I n 

oscillates (for an edge state) between and 1 , with every increment in i. However, a nonzero acts as a convergence 
factor. If it is quite small, but not zero, the ratio then goes smoothly to 1/2. Consequently, as Fig. 20 shows, 
'^'lo/^lo ^ 1/4 in the limit of small, but nonzero k. 

The effects illustrated in Fig. |20]- the three effects explicitly governed by k - are associated with the coupling 
between P and Q generated by T. Because this operator connects states with An = ±1, there is no large energy scale 
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TABLE III: Spectral property variations in H^^'^E) over 10 MeV 



Term 


Parameter 


1st Moment Shift (MeV) 


RMS Level Variation (MeV) 


Wave Function Overlaps 


{a\m 


K 


2.554 


1.107 


95.75-99.74% 


(a\V\l3) 


K 


0.272 


0.901 


99.35-99.82% 


{a\0\l3} 


K 


-0.239 


0.957 


99.51-99.99% 


{a\d{Em 


implicit 


0.135 


0.107 


99.95-100% 



associated with excitations. As the effects are encoded into a subset of the matrix elements, the overall scale of the 
K dependence on spectral properties is, at this point, still not obvious. 

This leaves us with one remaining term that, qualitatively, seems quite different, 

^ ^_^Q^ Q^^{aLo(|i^|), aNLoim, aNNLomi a^sLoim,-}- (42) 

Here the energy dependence is implicit, encoded in the parameters fitted to the lowest energy matrix elements of 
^e//^ The underlying potentials are dominated by strong, short-ranged potentials, much larger than nuclear binding 
energies. Thus the implicit ratio governing this energy dependence - \E\ vs. the strength of the hard-core potential - 
is a small parameter. For this reason one anticipates that the resulting energy dependence might be gentler than in 
the cases just explored. 

After repeating the fitting procedure over a range of energies, one obtains the results shown in Fig. |2T] Because 
the energy variation is quite small, results are provided only for the channels that contribute in low order, ^5*0, "^5*1, 
^Pi and ^Pj. The variation is very modest and regular, varying inversely with \E\ and well fit by the assumption 
(motivated by the form of V{E - QH)-^QV) 

ailQMeV) 

The variation is typically at the level of a few percent, over 20 MeV. The progression in the slopes within each channel, 
order by order, correspond to expectation: the lowest order terms, which account for the hardest part of the scattering 
in Q, have the weakest dependence on \E\. Comparisons between channels also reflect expectations. In the earlier 
discussion of naturalness, the rapid convergence in the ^5*0 channel, order by order, was consistent with very short 
range interactions in Q. Accordingly, a}j^Q varies by just 0.72% over a 10 MeV interval, and aj^j^Q by 1.10%. This 
channel contrasts with the ^Si channel, where convergence in the contact-gradient expansion is slower, consistent 
with somewhat longer range interactions in Q. For the "^S*! case one finds 2.64% variations in aLoi^Si) and 5.17% 
variations in aNLo^^Si) per 10 MeV interval. 

Are such variations of any numerical significance, compared to the explicit variations isolated in k? That is, if one 
were to determine a HOBET interaction directly from bound-state properties of light nuclei, would the neglect of 
this implicit energy dependence lead to significant errors in binding energies? One can envision doing such a fit over 
bound-state data spanning ~ 10 MeV, finding the couplings as a function of \E\, so that the error induced by using 
average energy-independent couplings aLo{\E\) can be assessed. These errors would reflect variations in the matrix 
elements to which these couplings are fit, following the procedures previously described. Such a study showed that 
only two channels exhibited drifts A in excess of 15 keV over 10 MeV, 

a^f? : A - ±21 keV af(^ : A ~ ±148 keV a^^^o : A - ±32 keV 

One concludes that the '^Si channel is, by a large factor, the dominant source of implicit energy dependence in the 
HOBET interaction. 

This allows one to do a more quantitative calculation that focuses on the most difficult channel {'^Si) and compares 
the relative sizes of the K-dependent and implicit energy dependences, as reflected in changes in the matrix i7'=-^^(|i;|). 
Thus this matrix is constructed at \E\ = 10 MeV and at \E\ ^ MeV (including the coupling to '^£'1), and changes in 
global quantities of that matrix over 10 MeV are examined: shifts in the first moment (the average eigenvalue) , the 
RMS shifts of levels relative to the first moment (related to the stability of level splittings), and eigenvalue overlaps. 
The four energy-dependent effects discussed here are separately turned on and off. Thus this exercise should provide 
a good test of the relative importance of these effects. The results are shown in Table [Hi] 

Despite the selection of the worst channel, ^Si, the implicit energy dependence is small, intrinsically and in compar- 
ison with the implicit energy dependence embedded in k. The implicit dependence in the first moment - a quantity 
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Energy (MeV) Energy (MeV) 

FIG. 21: The calculated energy dependence of derived coefficients for the contact-gradient expansion are indicated by the 
markers, for the various S — S and P — P channels. Over a 10 MeV interval typical of bound-state nuclear spectra, variations 
are typically at the few percent level. The continuous lines represent simple linear fits, a(10 MeV)/a(i5) = 1-1- a\E\, to the 
results. The fit is generally excellent. 



important to absolute binding energies - is 5% that of the explicit dependence in {a\T\(3). The RMS shifts in levels 
relative the first moment are at the ~ 1 MeV level for each of the implicit terms, but ^ 100 keV for the implicit term. 
Eigenfunction overlaps show almost no dependence on the implicit term, exceeding 99.95% in all cases: variations 
10-100 times larger arise from the analytical terms in k. 

Thus a simple representation of the HOBET effective interaction exists: 

• The requirements for a state of fixed \E\ are a series of short-range coefficients and a single parameter k that 
governs long-range corrections residing in Q, including certain terms that couple QV and QT. By various 
measures explored here, an N"^LO expansion is accurate to about a few keV 

• The K dependence found for a state of definite energy \E\ also captures almost all of the energy dependence 
resulting from varying \E\, the state-dependence in BH. Even in the most troublesome channel, calculations 
show that ~ 95% of the energy dependence associated with changes in \E\ is explicit. It appears that neglect 
of the implicit energy dependence would induce errors of < 100 keV, for a spectrum spanning ± 10 MeV. This 
kind of error would be within the uncertainties of the best ah initio excited-state methods for light (p-shell) 
nuclei, such as Green's function Monte Carlo 1 or large-basis no-core SM diagonalizations ilS]. 

• If better results are desired, the program described here can be extended to include the implicit energy depen- 
dence. The expansion around an average energy Eq 



V ^- QV = V 

E-QH^ 



Eo — QH Eq — QH Eq — QH 



QV 



generates the correction linear in E that is seen numerically. This second term is clearly quite small, explicitly 
suppressed by the ratio of scales discussed above. But, in any troublesome channel, the second term could be 
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represented by contact-gradient operators of low order, with the contribution suppressed by an overall factor of 

V. DISCUSSION AND SUMMARY 

One of the important motivations for trying to formulate an effective theory for nuclei in a harmonic oscillator 
basis is the prospect of incorporating into the approach some of the impressive numerical technology of approaches 
like the SM. Numerical techniques could be used to solve significant P-space problems, formulated in spaces, such as 
completed NTilo bases, that preserve the problem's translational invariance. 

My collaborators and I made an initial effort to construct a HOBET some years ago, using a contact-gradient 
expansion modeled after EFT. We performed a shell-by-shell integration, in the spirit of a discrete renormalization 
group method, but encountered several abnormalities connected with the running of the coefficients of the expansion. 
Subsequent numerical work in which we studied individual matrix elements revealed the problems illustrated in Fig. 
[3] These problems - the difficulty of representing Q-space contributions that are both long range and short range - 
are not only important for HOBET, but also are responsible for the lack of convergence of perturbative expansions of 
the effective interaction. Fig. [T| provides one example. In other work [5] we have shown that convergent expansions 
in the bare interaction (deuteron) or g-matrix (^He/^H) for H'^^^ do exist, if the long-range part of this problem is 
first solved, as we have done here. 

Thus the current paper returns to the problem of constructing a contact-gradient expansion for the effective in- 
teraction, taking into account what has been learned since the first, less successful effort. This paper introduces a 
form for that expansion that eliminates operator mixing, simplifying the fitting of coefficients and guaranteeing that 
coefficients determined in a given order remain fixed when higher-order terms are added. Thus the N^LO results 
presented here contain the results of all lower orders. The expansion is one in nodal quantum numbers, and is directly 
connected with traditional Talmi integral expansions, generalized for nonlocal interactions. 

Convergence does vary from channel to channel, but in each channel the order-by-order convergence is very regular. 
Each new order brings down the scale A at which deviations appear, and in each new order the Lepage plot steepens, 
showing that the omitted physics does have the expected dependence on higher-order polynomials in The 
channel-by-channel variations in convergence refiect similar behavior seen in EFT approaches, where the need for 
alternative power counting schemes has been noted to account for this behavior. From a practical standpoint, however, 
the N^LO results are effectively exact: in the most important difficult channel, ^Si, measures of the quality of the 
matrix W^^^ yielded results on the order of (1-3) keV. 

The summation done over QT yields a simple result, but still one that is quite remarkable in that long-range physics 
is governed by a single parameter k, that depends on the ratio of |£'| and fiLO. Despite all of the attractive analytic 
properties of the HO as a basis for bound states, its unphysical binding at large r has been viewed as a shortcoming. 
But the ladder properties of the HO in fact allow an exact summation of QT. It seems unlikely that any other 
bound-state basis would allow the coupling of P and Q by T to be exactly removed. That is, the HO basis may be 
the only one that allows the long-range physics in Q to be fully isolated, and thus subtracted systematically. In this 
sense it may be the optimal basis for correctly describing the asymptotic behavior of the wave function. Note, in 
particular, that the right answer is not going to result from using "improved" single-particle bases: k depends of \E\, 
not on single-particle energies of some mean field. 

The effects associated with k, are large, typically shifting edge-state matrix elements by several MeV, and altering 
spectral measures, like the first energy moment of H'^^^, by similar amounts. This dependence, if not isolated, destroys 
the systematic order- by-order improvement important to HOBET, as Fig. |3] clearly illustrates. 

The explicit energy dependence captured in n accounts for almost all of the energy dependence of H*^^^ {\E\). In 
more complicated calculations this dependence, in the BH formulation used here, generates the state-dependence that 
allows ET wave functions to have the proper relationship to the exact wave functions, namely that the former are 
the P-space restrictions of the latter. While in principle additional energy dependence important to this evolution 
resides in V{E — QH)~^QV and thus in the coefficients of the contact-gradient expansion, in practice this residual 
implicit energy dependence was found to be very weak. This dependence was examined channel by channel, and 
its impact on global properties of H'^^f{\E\) was determined for the most troublesome channel, '^S'l. Even in this 
channel, the impact of the remaining implicit energy dependence on H'^^^ {\E\) spectral properties such as the first 
moment, eigenvalue spacing, and eigenvalue overlaps, was found to be quite small compared to the explicit dependence 
isolated in k. This is physically very reasonable: QT generates nearest-shell couplings between P and Q, so that 
excitation scales are comparable to typical nuclear binding energies. Thus this physics, extracted and expressed as a 
function of k, should be sensitive to binding energies. In contrast, V{E — QH)~^QV involves large scales associated 
with the hard core, and thus should be relatively insensitive to variations in \E\. In the ^Si channel, the explicit 
dependence captured in k is about 20 times larger than the implicit energy dependence buried in the contact-gradient 
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coefficients. Numerically, the latter could cause drifts on the order of 100 keV over 10 MeV intervals. Thus, to an 
excellent approximation, one could treat these coefficients as constants in fitting the properties of low-lying spectra. 
Alternatively, the HOBET procedure for accounting for this implicit energy dependence has been described, and could 
be used in any troublesome channel. 

The weakness of the implicit energy dependence will certainly simplify future HOBET efforts to determine 
H'^f^{E, b, Ap) directly from data (rather than from an NN potential like vis)- Indeed, such an effort will be the next 
step in the program. The approach outlined here is an attractive starting point, as it can be shown that the states 
|5(k)) become asymptotic plane-wave states, when E is positive. Thus the formalism relates bound and continuum 
states through a common set of strong-interaction coefficients operating in a finite orthogonal basis. 

The relationship between current work and some more traditional treatments of the H^^-^ for model-based ap- 
proaches, like the SM, should be mentioned. Efforts like those of Kuo and Brown are often based on the division 
H = Ho + {V — Vo), where Hq is the HO Hamiltonian [T3]. Such a division would allow the same BH reorganization 
done here: QT and Q{Hq — Vq) are clearly equivalent. But in practice terms are, instead, organized in perturba- 
tion theory according to Hq, i.e., so that Green's functions involve single-particle energies. This would co-mingle 
the long- and short-range physics is a very complicated way. In addition, often the definitions of Q and P used in 
numerical calculations are not those of the HO: instead, a plane-wave momentum cut is often used, which simplifies 
the calculations but introduces uncontrolled errors. Either this approximation (plane waves are diagonal in T) or the 
use of perturbation theory (because of the co-mingling) would appear to make it impossible to separate long- and 
short-range physics correctly, as has been done here. 

Another example is Viow fc, in which a softer two-nucleon potential is derived by integration over high-momentum 
states [14]. This is a simpler description of Q than arises in bound-state bases problems, like those considered here: the 
division between P and Q is a specified momentum, and T is diagonal. There would be no analog of the K-dependence 
found in HOBET. However, HOBET and Viow k may have an interesting relationship. Effective operators for HOBET 
and for EFT approaches (which also employ a momentum cutoff) agree in lowest contributing order. When there are 
differences in higher order, it would seem that these differences must vanish by taking the appropriate limit, namely 
the limit of the HOBET (5(6, Ap) where 6 ^ oo while Ap/5 is kept fixed. This keeps the average (p^) of the last 
included shell fixed, while forcing the number of shells to infinity and the shell splitting to zero. Numerically it would 
be sufficient to approach this limit, so that Q resembles the plane- wave limit over a distance characteristic of the 
nuclear size. It is a reasonable conjecture that Viow k would emerge from such a limit of the HOBET H''ff{b,Ap). 
It would follow that all n dependence should vanish in that limit. It would be interesting to try to verify these 
conjectures in future work, and to study the evolution of the HOBET effective interaction coefficients as this limit is 
taken 

The state-dependence of effective interactions is sometimes treated in nuclear physics by the method of Lee and 
Suzuki. One form of Lee-Suzuki produces a Hermitian energy-independent interaction. While it is always possible 
to find such an H to reproduce eigenvalues, it is clear that basic wave function requirements of an effective theory ~ 
that the included wave functions correspond to restrictions of the true wave functions to P - are not consistent with 
such an H. 

Another form produces an energy-independent but nonHermitian H. This can be done consistently in an effective 
theory. However, the results presented here make it difficult to motivate such a transformation. It appears that the 
state-dependence is almost entirely attributable to the interplay between QT and QV, removed here analytically in 
terms of a function of one parameter k, which relates the bound-state momentum scale (in huj) to a state's energy. 
There is no obvious benefit in obscuring this simple dependence in a numerical transformation of the potential, given 
that the Lee-Suzuki method is not easy to implement. The physics is far more transparent in the BH formulation, and 
the self-consistency required in BH makes the use of an energy-dependent potential as easy as an energy-independent 
one. More to the point, the necessary k dependence is already encoded in the potential for a single state of definite 
energy \E\ - thus no additional complexity is posed by the state-dependence of the potential. 

I thank Tom Luu and Martin Savage for helpful discussions. This work was supported in part by the Office of 
Nuclear Physics and SciDAC, US Department of Energy. 

VI. REFERENCES 



[1] S. C. Pieper and R. B. Wiringa, Ann.Rev.Nucl.Part.Sci. 51 (2001) 53. 

[2] S. Weinberg, Phys. Lett. B251, 288 (1990), Nucl. Phys. B363, 3 (1991), and Phys. Lett. 295, 114 (1992). 



37 



[3] D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Lett. B 424 (1998) 390; S. R. Beane, P. F. Bedaque, W. C. Haxton, 
D. R. Philips, and M. J. Savage, in At the Frontier of Particle Physics, Vol. 1 (World Scientific, Singapore, 2001), p. 133; 
P. Bedaque and U. Van Kolck, Ann. Rev. Nucl. Part. Sci. 52 (2002) 334/ 

[4] R. B. Wiringa, V. Stoks, and R. Schiavilla, Phys. Rev. C51 (1995) 38. 

[5] W. C. Haxton and C.-L. Song, Phys. Rev. Lett. 84 (2000) 5484; C. L. Song, Univ. Washington Ph.D. thesis (1998) 

[6] W. C. Haxton and T. C. Luu, Nucl. Phys. A690, 15 (2001) and Phys. Rev. Lett. 89, 182503 (2002); T.C. Luu, S. Bogner, 

W. C. Haxton, and P. Navratil, Phys. Rev. C70, 014316 (2004). 
[7] S. Y. Lee and K. Suzuki, Phys. Lett. B91 (1980) 173; K. Suzuki and S. Y. Lee, Prog. Theor. Phys. 64 (1980) 2091. 
[8] L Talmi, Helv. Phys. Acta 25 (1952) 185. 
[9] R. Haydock, J. Phys. A7 (1974) 2120. 
[10] T. Luu, Ph. D. thesis. University of Washington (2003). 

[11] P. Lepage, How to Renormalize the Schrodinger Equation, lectures given at the VIH Jorge Andre Swieca Summer School 

(Brazil, 1997), nucl-th/9706029/. 
[12] P. Navratil and B. R. Barrett, Phys. Rev. C57 (1998) 562 andC59 (1999) 1906; P. Navratil, J. P. Vary, B. R. Barrett, 

Phys. Rev. Lett. 84, 5728 (2000). 
[13] T. T. S. Kuo and G. E. Brown, Nucl. Phys. A114 (1968) 241. 

[14] A. Schwenk, G. E. Brown, and B. Friman, Nucl. Phys. A703, 745 (2002); A. Schwenk, J. Phys. G31, S1273 (2005). 

[15] M. Valheres, T. K. Das, and H. T. Coelho, Nucl. Phys. A257 (1976) 389; also see G. Chen, Phys. Lett. A 328 (2004) 123. 



VII. APPENDIX: EFFECTIVE INTERACTION MATRIX ELEMENTS 



This appendix provides details on the evaluation of the modified HO states 

^ \nl). 



E-QT 



(A-1) 



for nodal quantum numbers n — 1,2,3,..., and for corresponding contact-gradient effective interaction matrix 
elements between such states. Closed- form expressions allow these matrix elements to be evaluated quickly to any 
order. Here two alternative evaluations are provided, one based on a harmonic oscillator expansion and one on the 
free Green's function. 



Harmonic Oscillator Expansion: The harmonic oscillator Green's function expansion is 
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where the gi are determined by a set of continued fractions generated from the ladder properties of the operator QT. 
In practice the sum can be truncated: numerical convergence is discussed below, in comparison with the Green's 
function approach. While this paper focuses on the simple example of the deuteron, the approach is more general : 
the relationship between the relative kinetic energy and the three-dimensional harmonic oscillator can be extended 
to the n-dimensional harmonic oscillator, with the hyperspherical harmonics replacing the spherical harmonics as 
eigenfunctions of the kinetic energy. The corresponding ladder properties for the harmonic oscillator in hyperspherical 
coordinates are known [15 : this is the essential requirement for the expansion. 

As discussed in the text, it is convenient to modify the usual contact-gradient expansion, to remove operator mixing 
and create an expansion in nodal quantum numbers. Each term O in the usual contact-gradient expansion is replaced 

by 



O^O 



(A-3) 



where r is the dimensionless Jacobi coordinate \ri — r2\/{bV^)- Gradients appearing in O are also defined in terms of 
this dimensionless coordinate. 

In each partial wave the lowest contributing operators are based on gradients, maximally coupled, acting on wave 
functions, with the result evaluated at r = 0. For example, for S, P, and D states 
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As the gradients are maximally coupled, all coupling schemes are equivalent. If one defines by (^'^)qQ the expressions 
with q gradients maximally coupled, the results of Eq. ( A-4 1 are examples of the more general formula 
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A form of this equation that will be used below is 
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Contact-gradient operators beyond the lowest contributing order involve acting on harmonic oscillator wave 
functions. One can quickly verify 
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so that 
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Thus by first using Eq. (A-8l and then applying Eq. (A-6 1, one finds the general expression for a contact-gradient 
operator of arbitrary order acting on a harmonic oscillator state 
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Equation (A-9 1 defines the needed matrix elements, evaluated below for each partial- wave channel contributing 
through N^LOfNote that if one wanted to write a potential, as opposed to partial-wave matrix elements of that 
potential that are given here, suitable projection operators could be inserted as needed. For example, the I — triplet 
and singlet channels could be distinguished by introducing the projection operators P{^Si) = (3 -I- cti • (72)/4 = S'^/2 
and (1 — (Ti • (T2)/4 = 1 — S^/2, respectively; the three triplet I = 1 channels could be distinguished from the singlet 
I — 1 channel and from each other by the projectors 
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and so on.] The matrix elements of Tabled] which all are scalar products of spin-spatial tensor operators, are of two 
types. One is diagonal in I, where O = O i ■ (5(f^ O; , with O; and O; spatial tensors. In this case 
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The second type of operator, O = t)i ■ S{f) [d^^2 ® [^(1) '^{'^)]2]i, enters for the off-diagonal triplet S — D,P — F, 
and D — G cases. 
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As both of these expressions reduce the matrix element to a product of terms like Eq. (A-9l, the needed N^LO 
matrix elements follow directly. For the '^Si ^ "^Si or ^ ^Sq channels, 
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For the ^ or ^Dj <-> ■^_Dj channels, recalling ~l5 j^j = ® '^)2M, 



{n'{l = 2S)J= IMTMt I e'''/^ 



y/2 



\n(l = 2S)J = IMTMt) = — ^ 

1 rz —z 



32 2 

15 71^ 



r[n' + 5/2]r[n + 5/2] 



-,1/2 



(7i'-l)!(n-l)! 



- 4[(n' - 1) + (n - l)]a]^^io 



(A-15) 



For the ^Ds ^ ^Gs channel, recalling = {{^ « '^)2 <S) (V ^ '^)2)4 «> (ct(1) <S) cr(2))i 
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For the ^ ^^^2 channel, recalling = ("^ (g) ("^ (g) '^)2)3 <g (a-(l) cr(2))5 
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And finally, for the ^F^ ^^F^ or ^Fj <-> ^Fj channels, recalling ~P\.i = <E) (V ig)'^)2)3M, 
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In each case the general matrix element of the effective interaction, for edge amd nonedge states, can then be 
expanded in terms of Eqs. ( A-13|[A^49 1, 
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Free Greenes Function Results: An equivalent set of results can be derived by making use of Eq. (36) 
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where the source term in the Green's function is 
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The inclusion of the last term makes this equation valid for non-edge as well as edge states: P\n + 1,^] = 1 if and 
only if |n + 1 Z) belongs to P. The reproduction of simple HO nonedge states is a helpful numerical check. Here 

The first task is to evaluate Ye^ on this wave function. The inclusion of e*" - which eliminates mixing 
among HO states and simplifies other HO expressions - is a bit of an annoyance in the Green's function case, generating 
a series of surface terms. One finds the generic result, analogous to Eq. (A-7|, 
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where each f{K^,r'^) is a polynomial that can be evaluated using standard gradient formul as for spherical harmonics. 
At N^LO /^4, is the highest contributing surface term. This form allows one to use Eq. ( A-7 1 to evaluate repeated 

operations of ^ . One can show, by expanding this expression around r = 0, that the leading order terms, which come 
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from the second line above and from the surface terms in Hne four and the following lines, are proportional to the solid 
harmonics r'^Yimi^r) , with corrections involving additional powers of r^. As the lowest order contributing operators, 
(^'^)qO, can annihilate only the former at the origin, in fact this expression effectively simplifies for contact-gradient 
purposes 
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For the cases of interest here (through N'^LO), the required nonzero polynomials are 
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105 + 105^2 + 21^4 + K^+ 4:1(1 + 1)(13 + 3^2) + 2Z(45 + 30^^ + 3k'^ + 41(1 + 1)) p=3 



fr\n\Q)+2lfl7\n^Q) = fJ,'(K^0)+2lfj:'(n\0) = fj,' (^^ 0) + 21 fJ/^K^ , 0) = 1 



fp=\^^',o)+2ifj\^^o) - 7+K^+2i rj:'(^^o)+2ifj:%K^o) = n+K^+21 



P=3/ 2 



fP=^(K\ 0) + 2lfj%K\ 0) = 57+18k"+k^+4/(/+1)+4?(6+2k^) 



The analog of Eq. (A-9l then becomes 



(^y(^\, e^'''K,(r)Yi^,(^r)\,^^ = -5i, ^ 



47r(2/ + l)!! 



f,(n^Q) + 2lfU>^\Q) 
fl(K\Q) + 2lfl,(n',Q) 



/+1/2 



( x^dx^Ki+if2(fix)RZi(x)^ + 2'+i 



2r[n + Z+ 1/2] 
7r(n- 1)! 



(^K^ + 3n + l- 3/2 + n, l)^yn(n + I + 1/2) + (n + I + l/2)P[n + 1,1]^ 



+16 



/^„(^^^0) + 2^/^3,(^^^0) 



r^,,(n\0) + 2lfl,,(n',0) 



(^(n - 1) [k^ + 3n + / - 5/2 + gi(~K^;n, I) ^/n(n + I + 1/2) ] 
+n(n + l+l/2)P[n+l,l] 



(n - l)(n - 2) [k^ + 3n + Z - 7/2 + gi(-K^;n, l)^Jn(n + I + 1/2) ] 
+n(n- l)(n + Z+ l/2)P[n + 1,/] 
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where 



(2Z + 1)!! 



r=0 



47r(2/+ 1)!! 



(A-26) 



has been used, and where the remaining integral can be evaluated using 
rw poo 1 

- / x^dx—Ki+i/2{Kx)R^^i{x) = -^2{n-iy.T[n + l + l/2] V 

Jo V _P 



(-2)^ 



^ m\(n - m)W\l + 3/2 + ml 

771—0 ^ / L / J 

(k^ + 3n + Z - m - 3/2 + gi{-K^\ nl)yjn{n + l + l/2){n - m) + n{n + I + l/2)P[n + 1, 1] 
' {V2Ky-' {l + iy. 



E 



i\{i-iy. 



-V2KT[m + 3/2+{l-i) /2] iFi[m + 3/2 + {I - i) /2; 3/2; /2] 
+T[m +! + (/- z)/2] iFi[m + 1 + (Z - i)/2; 1/2; k'^/2] . 



(A-27) 



With the effects of contact-gradient operators on wave functions of the form of Eq. (A-ll thus determined, these 
results can be plugged into Eqs. (A-ll I and (A- 12 1 to produce the needed expressions for matrix elements. 

The contact-gradient matrix elements needed through N'^LO have been evaluated using both the HO expansion and 
the free Green's function. The resulting agreement is a nice check. The sum over HO excitations is truncated at some 
N: at N'^LO the choice A^=400 will give results accurate to at least 0.01% for the most sensitive N'^LO operators, 
which depend on higher derivatives at the origin. The Mathematica script for this summation is simple and efficient, 
so there is no practical limit to the A^s that can be handled. Similarly, the Green's function expressions can be 
evaluated very easily: this is the recommended scheme for evaluating edge-state matrix elements for contact-gradient 
operat ors. 

Eq. ( A-21 1 is also an efficient way to generate the wave function at all values of r, which is needed for evaluations 
of edge-state matrix elements of the bare V. One finds 



Rniir) = y/2{n-iy.T[n + l + l/2] ^ 



(-1)" 

rn!(n-m)!r[/ + 3/2 + m] 



(n - m) + 3n - m + ; - 3/2 + gi{-K'^; nl)^yn{n + I + 1/2)^ +n{n + l+ l/2))P[n + 1, /] 
-^Ii+i/2{i^r)Gi[K,l,m,r] + -^Ki^i/2{Kr)G2[K,l,m,r] 
where 

^ ' ' ^ ^i\{l~iy (2k)*+i/2 ^ (2m+ 1 + ? - i - j)!jV72 

i — j — 

-V2{fi + r)T[l + j/2] iFi[l + j/2; 3/2; + rf I2\ + r[(l + j)/2] + j)/2; 1/2; + r) 2/2] 

{l-iy. 2™+('-*)/2 
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and 



G'2[K,Z,TO,r] = ^ V 



((-1)* + (-1)')k%/2 T[m + 3/2 + {I - i)/2] iFi [m + 3/2 + {I ^ i)/2; 3/2; k'^/2] 
+((-1)' - (-l)')r[m + i)/2] iFi [m + 1 + (/ - i)/2; 1/2; n^/2] 

2m+l+l-i 



E 



(2m+l + /-z)! r 2ra+i+i-i-j 
(2m+l + ;-i-j)!j! V2 



(-1)V«V2 {k-t) T[l+j/2] iFi[l+j/2;3/2;(«;-r)V2] 
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+(-l)'e-'='-\/2 + r) r[l + i/2] + i/2; 3/2; (k + r) V2] 

+{-lfe^^ r[(l + j)/2] iFi [(1 + j)/2; 1/2; - rf /2] 

-(-l)'e— r[(l + j)/2] + i)/2; 1/2; + rf /2] 



(A-30) 



